Semileptonic B and B s decays into orbitally excited charmed mesons 
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The BaBar Collaboration has recently reported products of branching fractions that include B 
meson semileptonic decays into final states with charged and neutral Di(2420) and 1)2(2460), two 
narrow orbitally excited charmed mesons. We evaluate these branching fractions, together with 
those concerning D ( *(2400) and Di(2430) mesons, within the framework of a constituent quark 
model. The calculation is performed in two steps, one of which involves a semileptonic decay and 
the other is mediated by a strong process. Our results are in agreement with the experimental 
data. We also extend the study to semileptonic decays of B„ into orbitally excited charmed-strange 
mesons, providing predictions to the possible measurements to be carried out at LHC. 
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I. INTRODUCTION 

Different collaborations have recently reported 
semileptonic B decays into orbitally excited charmed 
mesons providing detailed results of branching fractions. 
The theoretical analysis of these data, which include 
both weak and strong decays, offers the possibility for a 
stringent test of meson models. 

Moreover, an accurate determination of the \V c b\ and 
\V u b\ Cabbibo - Kobayashi - Maskawa matrix elements de- 
mands a detailed knowledge of semileptonic decays of b - 
hadrons. Decays including orbitally excited charmed me- 
son in the final state provide a substantial contribution 
to the total semileptonic decay width. Furthermore, a 
better understanding of these processes is also necessary 
in the analysis of signals and backgrounds of inclusive 
and exclusive measurements of 6-hadron decays. 

The Belle Collaboration |l|, using a full reconstruc- 
tion tagging method to suppress the large combinatorial 
background, reported data on the product of branching 
fractions B(B+ -> D**l + vi)B(D** -> D^tt), where, in 
the usual notation, I stands for a light e or /i lcpton, the 
Dq, D' Xl D\ and D 2 mesons are denoted generically as 
D**, and the D* and D mesons as DW. 

D** decays are reconstructed in the decay chains 
D** -> Z?*7T ± and D** -s- Dn ± . In particular, the D* Q 
meson decays only through the Dn channel, while the D[ 
and Di mesons decay only via D*ir. Both Dtt and D*ir 
channels are opened in the case of D\ . 

In the case of BaBar data @, Q the branching frac- 
tions B{D* 2 -> £)W-7r) include both the D* and D con- 
tributions. As they also provide the ratio B D / D (*) we 
estimate the D* and D contributions separately. The 
experimental results of both collaborations are given in 
Table H 

A similar analysis can be done in the strange sec- 
tor for the B s meson semileptonic decays. Here the 
intermediate states are the orbitally charmed-strange 
mesons, D** , and the available final channels are DK and 
D*K. The Particle Data Group (PDG) reports a value 



B(B° -> L> 5l (2536)-/i+^)B(D 8l (2536)- -> D*~K°) = 
2.4 ± 0.7 [|] based on their best value for B{b -> B°) 
and the experimental data for Bib — > B^.)B(B° S — > 
D sl (2536)~ /x+i/ /J )B(D s i(2536)- -> D*-R°) measured 
by the DO Collaboration [i|. 



All these magnitudes can be consistently calculated 
in the framework of constituent quark models because 
they can simultaneously account for the hadronic part of 
the weak process and the strong meson decays. In this 
context, meson strong decay has been described success- 
fully in phenomenological models, like the 3 Pn model Q 
or the flux-tube model [7J|, or in microscopic models 
(see Refs. @, Q). The difference between the two ap- 
proaches lies on the description of the qq creation vertex. 
While the 3 Pq model assumes that the qq pair is created 
from the vacuum with vacuum quantum numbers, in the 
microscopic model the qq pair is created from the in- 
terquark interactions already acting in the model. Both 
approaches will be used here to evaluate the strong de- 
cays. As for the weak process the matrix elements factor- 
izes into a leptonic and a hadronic part. It is the hadronic 
part that contains the nonperturbative strong interaction 
effects and we shall evaluate it within a constituent quark 
model (CQM). We will work within the CQM of Ref. 
which successfully describes hadron phenomenology and 
hadronic reactions |ldH13l | and has recently been applied 
to mesons containing heavy quarks in Refs. [141 Il5|. 



The paper is organized as follows: In Sec. [TT] we in- 
troduce the model we have used to get the masses and 
wave functions of the mesons involved in the reactions 
mentioned above. In Sees. IIIII and IIVI we study the 
semileptonic and strong decay mechanisms, which consti- 
tute the two steps of the processes under study. Finally, 
we present our results in Sec. [V] and give some conclusions 
in Sec. [VE 



Belle [1] (xlO" 3 ) BaBar [2, 3] (xlCT 



DS(2400) 

B(B+ -> Dl°l + vi)B{D* ° -> D-^ 
B(B° -» D*~l + i>i)B(D*- ->■ DV 



2.4 ±0.4 ±0.6 
2.0 ±0.7 ±0.5 



2.6 ±0.5 ±0.4 
4.4 ±0.8 ±0.6 



Di(2430) 






'o 



► D^l + v{)B{D 



D*--k+) 

>D*°TT-) 



<0.7 
<5 



2.7 ±0.4 ±0.5 
3.1 ±0.7 ±0.5 



Di(2420) 

B(B+ 
B(B° 



y Dll+vi)B{Dl -> D*-tt+) 
D-I+v^bId- -+D* tt ~) 



4.2 ±0.7 ±0.7 
5.4 ±1.9 ±0.9 



2.97 ±0.17 ±0.17 
2.78 ±0.24 ±0.25 



Da(2460) 

B(B H 
B(S H 



Dfl+ Vl )B(Df 
D?l+v{)B{D? 



B(B° 
B(B° 
B{B° 



DX~l + vi)B{D* 2 
, 2 >-vi)B{pl- 






d;-i+^)b(d;- 



D-TT+) 

D*-7v+) 



B 



d/dM 



2.2 ±0.3 ±0.4 
1.8 ±0.6 ±0.3 
4.0 ±0.7 ±0.5 

2.2 ±0.4 ±0.4 

<3 

< 5.2 

0.55 ± 0.03 



1.4±0.2±0.2 ( * ) 

0.9±0.2±0.2 W 

2.3 ±0.2 ±0.2 

1.1±0.2±0.1 ( * ) 

0.7±0.2±0.1 ( * ) 

1.8 ±0.3 ±0.1 

0.62 ± 0.03 ± 0.02 



TABLE I. Most recent experimental measurements reported by Belle and BaBar Collaborations. I stands for a light e or fj, 
lepton. The symbol (*) indicates results estimated from the original data by using B D , D (,). 



II. CONSTITUENT QUARK MODEL 



Spontaneous chiral symmetry breaking of the QCD La- 
grangian together with the perturbative one-gluon ex- 
change (OGE) and the nonperturbative confining inter- 
action are the main pieces of potential models. Using 
this idea, Vijande et al. [10| developed a model of the 
quark-quark interaction which is able to describe meson 
phenomenology from the light to the heavy quark sector. 
We briefly explain the model below. Further details can 
be found in Ref. [13]. 

One consequence of the spontaneous chiral symmetry 
breaking is that the nearly massless 'current' light quarks 
acquire a dynamical, momentum-dependent mass M(j>) 
with M(0) sw 300 MeV for the u and d quarks, namely, 
the constituent mass. To preserve chiral invariance of 
the QCD Lagrangian new interaction terms, given by 
Goldstone boson exchanges, should appear between con- 



stituent quarks. 

A simple Lagrangian invariant under chiral transfor- 
mations can be derived as [la] 



£ = VOt^-MC/ 75 )^, 



(1) 



where If 5 = exp(z7r a A a 75//,r), 7r a denotes the pseu- 
doscalar fields (7?, K, r]s) and M is the constituent quark 
mass. The momentum-dependent mass acts as a natural 
cutoff of the theory. The chiral quark-quark interaction 
can be written as 

V n (fij) = Vg (nj) + V q T q (fy) + V q f (f tJ ) , (2) 

where C, T and SO stand for central, tensor and spin- 
orbit potentials. The central part presents four different 
contributions, 

V q ° q (nj ) = K c (fa ) + V? (n, ) + V° (r tJ ) + v v c (r l3 ) , (3) 

given by 



V?(fa) 
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A 3 
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(oW,)[cos0 p (Af-A|)- 



(4) 



where Y(ir) is the standard Yukawa function defined by 
Y(x) = e~ x J x. We consider the physical v. meson instead 
of the octet one and so we introduce the angle 9 p . The A a 
are the SU(3) flavor Gell-Mann matrices, m, is the quark 
mass and m^, ttik and m v are the masses of the SU(3) 
Goldstone bosons, taken from experimental values. m a 
is determined through the PCAC relation m 2 ~ m 2 + 
4m 2 d [17|. Finally, the chiral coupling constant, g c h, is 
determined from the irNN coupling constant through 



a 2 

9ch 

4tt 25 4tt 



9 91nn m u,d 



'N 
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which assumes that flavor SU(3) is an exact symmetry 
only broken by the different mass of the strange quark. 

There are three different contributions to the tensor 
potential 



V£(fij) = V?(? tj ) + Vg(? tj ) + V?(fij), 



(6) 



(5) given by 
I 



Vj (fj) 



A 



4-7T Ylrriimj A 2 — m 2 



VI [fj) 



A 2 



• f (p..) — 9ch mK 

An 12rriimj A\ — m\ 
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H{m K rij) - —^-H(A K r ij ) 

K 



Vfifj) 



9ch m2 n 



A 2 



Yk YlmiUij A 2 — m 



2"*V 



A 3 



H(mr,rij) ^H(A v r ij ) 



a=l 

7 

Sij / yJK ' a°)i 

a=4 

Sij [cos0 p (\f -A?) - sin 9 p ] 



(7) 



Sij = i{<Ji ■ T%j){<Jj ■ r%j) — <Ti ■ <Jj is the quark tensor 
operator and H(x) = (1 + 3/x + 3/x 2 )Y(x). 

Finally, the spin-orbit potential only presents a contri- 
bution coming from the scalar part of the interaction 



Y 



3 (f ;,. ) = K so (f )= 9l h ml A 2 

J J 4tt 2m l m j A 2 , - ml 


r A 3 1 
G(m <T r lJ ) |G(A CT r^) 


(L-S). 



(8) 



In the last equation G(x) is the function (l + l/x)Y(x)/x. 
Beyond the chiral symmetry breaking scale one expects 
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the dynamics to be governed by QCD perturbative ef- 
fects. In this way one-gluon fluctuations around the in- 
stanton vacuum are taken into account through the qqg 
coupling 



WAna^j^X ^, 



(9) 



with A c being the SU(3) color matrices and G% the gluon 
field. 

The different terms of the potential derived from the 
Lagrangian contain central, tensor, and spin-orbit con- 
tributions and are given by 



V§ GE (nj) = i a s (A l c • Ap 



* / OGE( r *j) — 



^oge(%) 



1 a s 

16 TOjTOj 
1 Q S 

16 m 2 m 2 
1 j 



6rriimj 



(di -d 3 y 



-rtj/robi) 



(Xt ■ A?) 
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e -ry/r g (|t) / ^ 



e -nj/r g (n) 
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ij 



^s(M) 



4 3r2(^) r^r^) 



((m, + m^) 2 + 2m l m : j)(S + ■ L) + (m 2 - mf)(S- ■ L) 



(10) 



where S± — \{3i ± 3j)- Besides, r (fi) = fo 1 ^- and 
r g (/^) = fg 1 ^- are regulators which depend on /uy, the 
reduced mass of the qq pair. The contact term of the 
central potential has been regularized as 



S(fi 



-ru/ro 



Airr 2 , 



(11) 



The wide energy range needed to provide a consistent 
description of light, strange and heavy mesons requires an 
effective scale-dependent strong coupling constant. We 
use the frozen coupling constant of Ref. [10( 



,M = 



Of) 



In 



(12) 



r 



in which fi is the reduced mass of the qq pair and ao, 
Ho and Ao are parameters of the model determined by a 
global fit to the meson spectra. 

Confinement is one of the crucial aspects of QCD. 
Color charges are confined inside hadrons. It is well 
known that multigluon exchanges produce an attractive 
linearly rising potential proportional to the distance be- 
tween quarks. This idea has been confirmed, but not rig- 
orously proved, by quenched lattice gauge Wilson loop 
calculations for heavy valence quark systems. However, 
sea quarks are also important ingredients of the strong 
interaction dynamics. When included in the lattice cal- 
culations they contribute to the screening of the rising 
potential at low momenta and eventually to the breaking 
of the quark-antiquark binding string. This fact, which 
has been observed in n/ = 2 lattice QCD [18| . has been 
taken into account in our model by including the terms 



coiwy) - 
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.•£)_ 


} 



2a s ) +4miTO_y(l - a s ))(S+ ■ L) 



(13) 



where a s controls the mixture between the scalar and 
vector Lorentz structures of the confinement. At short 
distances this potential presents a linear behavior with 
an effective confinement strength a = — a c /i c (A? • Xj) 
and becomes constant at large distances with a threshold 
defined by 



Khr = {-a c + A}(A?-A?). 



(14) 



No qq bound states can be found for energies higher 
than this threshold. The system suffers a transition 
from a color string configuration between two static color 
sources into a pair of static mesons due to the breaking of 
the color string and the most favored decay into hadrons. 

Among the different methods to solve the Schrodinger 
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equation in order to find the quark-antiquark bound 
states, we use the Gaussian Expansion Method 19] be- 
cause it provides enough accuracy and it makes the subse- 
quent evaluation of the decay amplitude matrix elements 
easier. 

This procedure provides the radial wave function solu- 
tion of the Schrodinger equation as an expansion in terms 
of basis functions 



R a (r) = J2 <€i{r), 



(15) 



where a refers to the channel quantum numbers. The 
coefficients, c", and the eigenvalue, E, are determined 



Quark masses 


m„ (MeV) 


313 




m s (MeV 


555 




m c (MeV 


1763 




m b (MeV) 


5110 


Goldstone Bosons 


m^ (fin -1 


0.70 




mc (fm - 


3.42 




7TCk (fm -1 


) 2.51 




m n (fm _1 


2.77 




A^ (fm" 1 ) 


4.20 




K (fm" 1 ) 


4.20 




A K (fm" 1 


4.21 




A„ (fm -1 ) 


5.20 




<? C \/47T 


0.54 




M°) 


-15 


OGE 


ao 


2.118 




Ao (fm" 1 ) 


0.113 




Mo (MeV) 


36.976 




r (fm) 


0.181 




r s (fm) 


0.259 


Confinement 


a c (MeV) 


507.4 




iu c (fm -1 ) 


0.576 




A (MeV) 


184.432 




a s 


0.81 



TABLE II. Quark model parameters. 



from the Rayleigh-Ritz variational principle 



E 



\-*-rin ^^n'nJ^n' 



yaa a _ q 



, (16) 



where T°, n , N", n and V%™ are the matrix elements of 
the kinetic energy, the normalization and the potential, 
respectively. T", n and N", n are diagonal whereas the 
mixing between different channels is given by V$" . 

Following Ref. [19J , we employ Gaussian trial functions 
with ranges in geometric progression. This enables the 
optimization of ranges employing a small number of free 
parameters. Moreover, the geometric progression is dense 
at short distances, so that it allows the description of 
the dynamics mediated by short range potentials. The 
fast damping of the gaussian tail is not a problem, since 
we can choose the maximal range much longer than the 
hadronic size. 

Table [III shows the model parameters fitted over all 
meson spectra and taken from Refs. fl(X[l4j]. 



excited charmed mesons. In the nonstrange sector, this 
has been studied before within heavy quark effective 
theory (HQET) in Refs. 0, HJ. There, only relative 
branching ratios could be predicted and their results 
depended on the approximation used and on two un- 
known functions, t\, T2, that describe corrections of or- 
der Aqcd/^q. Only the ratio T^»/Vd**, semileptonic 
decay rate with a helicity D** final meson over total 
semileptonic decay rate to that meson, seemed to be sta- 
ble in the different approximations. We shall comment 
on this below. 

In the context of nonrelativistic constituent quark 
models, the state of a meson is given by 



\M,XP) 



NR 



E 



( _l)l/2- 



d 3 p 
x <& A 2 } (P) I 9, «i Pi > Iff, "a Vi ) > (17) 



where P is the three-momentum of the meson and A is 
the spin projection in the meson center of mass. The 
vector p is the relative momentum of the qq pair, p\ = 



-P—pax\dp2 = 



V 2 



-P+p are the momenta of 



m fl +m f2 f r* m fl +m f2 

the antiquark and the quark, respectively, a.\ and «2 are 
the spin, flavor and color quantum numbers. (E(pi),pi) 
are the four-momenta and rrii are the quark masses. The 
factor (— l) 1 / 2 - 151 is included in order that the antiquark 
spin states have the correct relative phase. 

The normalization of the quark-antiquark states is 

(a'p'\ap) = 8 a ,, a (2tt) 3 2E f (p) 5{p' -0), (18) 

and the momentum space wave function (f>ax,ai \{p ) nor_ 
malization is given by 



d 3 P £(*&£W#££>(?) = * 



A', A- 



(19) 



Finally, the normalization of our meson states is 

l (M,\'P'\M 1 XP) NR = dy, x (27r) 3 d(P / - P). (20) 



NR\ 



In the decay we have a b — > c transition at the quark 
level and we need to evaluate the hadronic matrix ele- 
ments of the weak current 



III. WEAK DECAYS 



J£ c (0)= ^(0)7^-75)^(0). 



(21) 



In this section, we give an account of the semilep- 
tonic decays of the B {B or B s ) meson into orbitally 



The hadronic matrix elements can be parameterized in 
terms of form factors as 



(D(0+),\P D \J» c (0)\B(0-),P B ) = P^F + {q 2 )+q^F_{q 2 ), 



(D(l+),XP D \Jf(0)\B(0-),P B ) 



e^ aP eY* ) (P D )P a qPA(q 2 ) 



m B + mo 
- i I (m B - m D )e* x ^{P D )Vo(q 2 ) 
(D(2+),XP D \J" C (0)\B(0-)P B ) =e^ a tflUP D )P s P a q l3 T A {q 2 ) 



m B + mo 



i 



{^X)^(PD)P S Tx(q 2 ) + P"P S e* w „ s (P D ) [P»T 2 (q 2 ) + q»T 3 (q 2 )} } . 



(22) 



In the expressions above, P = P B + P B and q = P B — P B , 
P B and Pd being the meson four- moment a. m B and mo 
are the meson masses, <^ LVa P is the fully antisymmetric 
tensor, for which the convention e 0123 = +1 is taken, and 

e (x)/j,(P) an( i € (\)iiv{P) are the polarization vector and 
tensor of vector and tensor mesons, respectively. The 
meson states in the Lorentz decompositions of Eq. 
are normalized such that 



(M,X'P'\M,XP) 



J\'\ 



{27r) 3 2E M (P)6{P' - P). (23) 



where Em (P) is the energy of the M meson with three- 
momentum P. Note the factor 2Em difference with re- 
spect to Eq. (|2Ull . 

The form factors will be evaluated in the center of mass 

I 



r 



of the meson, taking q in the z direction, so that 
P B = and Pd = —q = — \q\k, with k representing the 
unit vector in the z direction. We have taken the phases 
of the states such that all form factors are real. F + , F_ , 
A, Vo, V+, V- and T\ are dimensionless, whereas T 2 , T 3 
and T4 have dimension of E~ 2 . Defining vector Vt(\q |) 
and axial A^(|<f|) matrix elements such that 



V£{\q\) = (M F ,X- \q\k\J b v c ^0)\Mi,6) , 



A»(\q*\) = (M F ,X- \q\k\J* A c »(0)\Mi,0) , 



(24) 



we have for a —¥ + decay, that the form factors are 
given in terms of vector and axial matrix elements as 



F + {q 2 



F-(q 2 
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(25) 



I 

In the case of a _ — > 1 + transition, the corresponding expressions for the form factors are 
I 
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Finally, the form factors for a — > 2 + transition are given by the relations 



Tl(q 2 ) = - 



. 2m D Al 



ATx=+i(\q 




,2mo( E D (-q)(E D (-q) - m B ) \ t \ 
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(27) 



The CQM evaluation of the vector and axial matrix ele- 
ments V^ (| q | ) and j4 a (| q | ) can be found in the Appendix. 
For a B meson at rest and neglecting the neutrino 
mass, we have the double differential decay width 



d 2 T 
dq 2 dxi 
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\V bc \ 2 W\? 
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8tt 3 2m B 

■H a p{P B ,P D )C aP {p hPy ), 



(28) 



where xi is the cosine of the angle between the final me- 
son momentum and the momentum of the final charged 
lepton measured in the lepton-neutrino center of mass 
frame. G F = 1.16637(1) x 1CT 5 GeV -2 is the Fermi 
constant jj, mi is the charged lepton mass, A(a, b, c) = 
(a + b — c) 2 — Aab and V& c is the be element of the Cab- 
bibo - Kobayashi - Maskawa matrix for which we shall use 
Vbc = 0.0413. H a p and C af3 represent the hadron and 
lepton tensors. P B , Pd, Pi and p v are the meson and 
lepton momenta. 

Working in the helicity formalism of Ref. [22| and after 
integration on xi we have 



dq 2 "Stt 3 ' bc 
x(H v 



M 2 
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2 ,2 



A 1 / 2 ^ 2 ,'' 



L Bi • 



ll) 



12m B q 2 

Hl + Hu 



Hr 



2m B 



(29) 



where the suffixes U, L, S stand for the unpolarized- 
transverse, longitudinal and scalar components of the 
hadronic tensor, and H — IpzH. Integrating over q 2 
we obtain the total decay width that can be written as 



r = r 



u 



•r. 



u 



r s , 



with Tj and Tj partial helicity widths defined as 
- 2 m 2 ) 2 X 1 / 2 {q 2 ,m 2 Bl m 2 D ) 



r -'*#^l 2i W 



2m B 



(30) 

Hj 

(31) 
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and similarly for Tj in terms of Hj. The evaluation 
of the different form factors, and thus of the different 
helicity amplitudes of the hadronic tensor, has been done 
following Ref. [H[. 



IV. STRONG DECAYS 

Meson strong decay is a complex nonperturbative pro- 
cess that has not yet been described from QCD first prin- 
ciples. Instead, several phenomcnological models have 
been developed to deal with this topic, the 3 Pq p, the 
flux-tube [3], and the Cornell [8|,[9j models being the most 
popular. 

Some models describe the decay process assuming that 
the extra quark-antiquark pair is created from the vac- 
uum. This is the case of the 3 Pq model, which borrows 
its name from the quantum numbers of the created pair, 
or the flux-tube model, which in addition to the creation 
vertex incorporates the overlaps between the color flux 
tubes of the initial and final states. 

To address a more fundamental description of the de- 
cay mechanism, one has to describe hadron strong decays 
in terms of quark and gluon degrees of freedom. However, 
there has been little previous work in this area. Two dif- 
ferent examples are the study of open-charm decays of cc 
resonances by Eichten et al. [8j, who assumed that the 
decays are due to pair production from the static part of 
a Lorentz vector confining interaction, and the study of a 
few strong decays in the light sector by Ackleh et al. [9| , 
where the qq pair production comes from the one-gluon 
exchange and a scalar confining interaction. 

As we mentioned in the introduction, we shall use both 
the 3 Pq model and a microscopic one, resembling those of 
Refs. Q and Q, that originates from the different inter- 
action pieces present in our interquark potential. These 
two approaches to meson production are introduced in 
the following subsections. 
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A. The S P model 

It was first proposed by Micu Q and further devel- 
oped by Le Yaouanc et al. [24]. To describe the meson 
decay process A — » B + C, the 3 Pq model assumes that a 
quark-antiquark pair is created with vacuum J PC = ++ 
quantum numbers. The created qq pair together with the 
qq pair present in the original meson regroup in the two 
outgoing mesons via a quark rearrangement process. 

The interaction Hamiltonian which describes the pro- 
duction process is given by [9J 



H, 



g / d 3 xip(x)ip(x) 



(32) 



where g is related to the dimensionless constant giving 
the strength of the qq pair creation from the vacuum 
as 7 = 2^T' m 9 being the mass of the created quark. 
Note that the operator g'tptp leads to the decay (qq)A — > 
(qq)b + (qq)c through the attf term. 



B. The microscopic model 

In microscopic decay models one attempts to describe 
hadron strong decays in terms of quark and gluon degrees 

I 



of freedom. The quark-gluon decay mechanism should 
give similar predictions to the reasonably accurate 3 Pq 
model and should determine the strength of the qq pair 
creation, 7, of the 3 Pq model in terms of more fundamen- 
tal parameters. 

Following Ref. [9J, the strong decays should be driven 
by the same interquark Hamiltonian which determines 
the spectrum, the one-gluon exchange, and the confin- 
ing interaction appearing as the kernels. These inter- 
actions and their associated decay amplitudes are un- 
doubtedly all present and should be added coherently. 
We already mentioned that our constituent quark model 
for the heavy quark sector has a one-gluon exchange 
term and a mixture of Lorentz scalar and vector confin- 
ing interactions. This completely defines our microscopic 
model for strong decays. Unlike previous works we use 
a screening confinement interaction and also a mixture 
between scalar and vector Lorentz structures, which is 
already fixed. 

The Hamiltonian of the interaction can be written as 



Hj = ~ J d 3 xd 3 y J a (x)K(\x - y\)J a (y). 



(33) 



The current J a in Eq. (|33|) is assumed to be a color octet. 
The currents, J, with the color dependence A a /2 factored 
out and the kernels, K(r), for the interactions are 



• Currents 



j(x) = ip(x) r i/j(x) 



■0(x)X^(x) Scalar Lorentz current, 

ijj(x) 7 ip(x) Static part of vector Lorentz current, 

V>(x) 7 V>(^0 Spatial part of vector Lorentz current, 



(34) 



• Kernels 



K{r) 



-4a s [-o c (l - e-»") + A] 



Confining interaction, 
Color Coulomb OGE, 
Transverse OGE. 



(35) 



For the Lorentz vector structure of the confinement 
we use K(r) = ±(1 — a s )4 [— a c (l - e _Alc '') + A], where 
± refers to static and transverse terms, respectively. We 
refer, following Ref. Q, to this general type of interac- 
tion as a JKJ decay model, and to the specific cases 
considered here as sKs, j°Kj° and j T Kj T interactions. 



C. Strong decay width 

The total width is the sum over the partial widths 
characterized by the quantum numbers Jbc and I 



^a^bc— 2^, ^a^bc(JbcJ) 



(36) 



Jbc, l 



The wave functions for the mesons involved in the reac- 
tions are the solutions of the Schrodinger equation using 
the Gaussian Expansion Method mentioned above. De- 
tails of the resulting matrix elements for different cases 
are given in Ref. [25j . 



where 



T A -+bc{Jbc, I) = 2tt f dk Q S(E A - E BC )\M A ^Bc(ko) 



37) 



and M.A^Bc(ko) is calculated according to Refs. 25, 2 



Using relativistic phase space, we arrive at 

r A ^Bc(JBC,l) = 2TT^^\M A ^Bc(ko)\ 2 , (38) 
m A k 



where 



fco = 



\ 1 / 2 {m 2 A ,m 2 B: m 2 c ) 
2m A 



(39) 



is the on shell relative momentum of mesons B and C. 



V. RESULTS 

For the low-lying positive parity excitations, any quark 
model predicts four states that in the 2S+1 Lj basis cor- 
respond to 1 Pi, 3 Pq, 3 Pi and 3 P2. As charge conjugation 
is not well defined in the heavy-light sector, l P\ and 3 Pi 
states can mix under the interaction. 

In the infinite heavy quark mass limit, heavy quark 
symmetry (HQS) predicts two degenerated P-wave me- 
son doublets, labeled by j q — 1/2 with J p = + ,l + 
(|l/2,0+),|l/2,l+)) and j q = 3/2 with J p = l+,2+ 
(|3/2, 1 + ), |3/2, 2 + )). In this limit, the meson properties 
are governed by the dynamics of the light quark, which is 
characterized by its total angular momentum j q = s q +L, 
where s q is the light quark spin and L the orbital angular 
momentum. The total angular momentum of the meson 
J is obtained coupling j q to the heavy quark spin, sq. 

Moreover, in the infinite heavy quark mass limit the 
strong decays of the Dj (j q = 3/2) proceed only through 
D-waves, while the Dj (j q = 1/2) decays happen only 
through S- waves [27[. The D-wave decay is suppressed 
by the barrier factor which behaves as q 2L+1 where q 
is the relative momentum of the two decaying mesons. 
Therefore, the states decaying through D- waves are ex- 
pected to be narrower than those decaying via iS-waves. 

A change of basis allows to express the above states 
in terms of the 2S+1 Lj basis, by recoupling angular mo- 
menta, as 



11/2,0+) = 


= +| 3 Po) 




|l/2,l + > = 


-+\fl\ lp ^ + \ 


/I** 


|3/2,l+) = 


-V!' 1 *^ 


fa* 


13/2,2+) = 


= + |'ft) 





(40) 



where in the 2 +1 Lj wave functions we couple heavy and 
light quark spins, in this order, to total spin S. 

In the actual calculation the ideal mixing in Eq. (j40|) 
between 1 Pi and 3 Pi states changes due to finite charm 
quark mass effects. Our CQM model predicts the mixed 
states shown in Table IIII1 which are very similar to the 
HQS states. This is expected since the c- quark is much 
heavier (m c = 1763 MeV) than the light (m n = 313 MeV) 
or strange (m s = 555 MeV) quarks. Note that now we 



3 Po 
'Pi 
3 A 

3 P 2 

1/2,0+ 
1/2,1+ 
3/2, 1+ 
3/2,2+ 



Wo 
+, 1.0000 



A 

-, 0.5903 
+, 0.4097 



A 



0.4097 
0.5903 



Dl 



+, 0.99993 



+, 1.0000 



+, 0.0063 -, 0.9937 
+, 0.9937 +, 0.0063 



0.99993 



Ao 

3 P +, 1.0000 

'Pi 

3 Pi 

3 P 2 

1/2,0+ +, 1.0000 
1/2,1+ 
3/2, 1+ 
3/2,2+ 



Ai 



Ai 



0.7210 -, 0.1880 
0.2770 -, 0.5570 



d:. 



0.99991 



0.0038 
0.9942 



0.7390 
0.0060 



+, 0.99991 



TABLE III. Probability distributions and their relative phases 
for the four states predicted by CQM in the two basis de- 
scribed in the text. In the 1+ strange sector the effects of 
non-gg components are included, see text for details. 



have mixing, even if small, between the 3 P2 and 3 P2 par- 
tial waves in 2+ mesons. This is due to the OGE tensor 
term. 

In Ref. [15| we have studied the J p = 1+ charmed- 
strange mesons, finding that the J p = 1+ P s i(2460) 
has an important non-qq contribution whereas the J p = 
1+ -D s i(2536) is almost a pure qq state. The presence 
of non-qq degrees of freedom in the J p = 1+ charmed- 
strange meson sector enhances the j q — 3/2 component 
of the _D s i(2536). This wave function explains most of 
the experimental data, as shown in Ref. [15], and it is 
the one we shall use here. For this sector only the qq 
probabilities are given in Table Mil 



A. B semileptonic decays into D** mesons 



Semileptonic B — >■ D* ) {24QQ)lvi decay 



The 
B(B+ 
B(B° 



measured branching 
D* l+ Vl )B(D* 



fractions 



D* -i+ n )B(D* - 



P°7T-). 



are 

~7r+) and 
The meson 



Pq(2400) has J p = 0+ quantum numbers and, therefore, 
due to parity conservation, it decays only into Dir, so that 
we have B(D^° -» p-7r+) = B{D*^ -> P 7r-) = 2/3 
coming from isospin symmetry. 

Table IIVI shows the different helicity contributions to 
the semileptonic width. In both cases the dominant con- 
tribution is given by Yl while the rest are negligible. The 
difference between the semileptonic width of the charged 
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B+ -y Dfl + m 


B° -> D^-l+n 


rv 


0.00 


0.00 


IV 


0.00 


0.00 


r L 


1.30 


1.16 


r L 


6.83 x 10" 7 


6.45 x 10 -7 


r s 


2.05 x 10~ 6 


1.93 x 10~ 6 


r 


1.30 


1.16 



TABLE IV. Helicity contributions and total decay widths, in 
units of 10 -15 GeV, for the Dq meson. 





B+ -> Dll+vt 


B° ->■ D^l + vi 


IV 


0.38 


0.38 


IV 


1.94 x 10" 8 


1.93 x 10~ 8 


IV 


1.17 


1.16 


r L 


7.16 x 10~ 7 


7.15 x 10~ 7 


r s 


2.17 x 10~ 6 


2.17 x 10~ 6 


r 


1.55 


1.54 



TABLE VI. Helicity contributions and total decay widths, in 
units of 10~ 15 GeV, for the Di meson. 





B+ ->■ D^l+vi 


B° -4 Dfl+n 


r v 


0.23 


0.23 


IV 


1.35 x 10~ 8 


1.35 x 10~ 8 


r L 


0.56 


0.56 


r L 


4.12 x 10 -7 


4.12 x 10~ 7 


r s 


1.27 x 10" 6 


1.27 x 10~ 6 


r 


0.79 


0.80 



TABLE V. Helicity contributions and total decay widths, in 
units of 10~ 15 GeV, for the D[ meson. 



Wc have in this case the product of branching fractions 

B{B+ -> D°l + vi)B{D^ -> D*-ir+) = 1.32 x 10~ 3 , 

B(B° -> Dfl+^BiDf -> D*°ir-) = 1.23 x 1(T 3 . 

(42) 

which are a rough factor of 2 smaller than the results 
from the BaBar Collaboration [2j . 

3. Semileptonic B -+ Di(2420)lui decay 



and neutral B meson is due to the large mass difference 
between the Dq and Dq mesons for which we take the 
masses reported in Ref. [4] . 

Figure Q] shows the q 2 dependence in the form factors 
and in the differential decay width for B(B + — > Dq°1 + vi), 
panels (a) and (b), respectively. Similar results (not 
shown) are obtained for the B(B° —t Dq~1 + vi) case. 

The final results for the product of branching fractions 



B(B+ -> Dl°l + vi)B{Dl° -> D-TT+) 
B(B° -> D*-l + Vi)B(D*- -> D\~) 



2.15 x 1(T 3 , 
1.80 x 10~ 3 , 



(41) 



which compare very well with Belle data |l[ . 



2. Semileptomc B — > D[ (2430) /i^ decay 

The only Okubo - Zweig - Iizuka (OZI)-allowed decay 
channel for the D[ meson is D[ — > D*tt so that 
isospin symmetry predicts a branching fraction B(D[ — > 
D*n ± ) = 2/3. 

Table |V] shows the different helicity contributions to 
the semileptonic width of B + —} D®l + vi and B° —t 
D^l + vi calculated in the framework of the CQM. In this 
case, Tu and IV are of the same order of magnitude and 
give the total semileptonic decay rate. 

Panels (a) and (b) of Fig. [2] show the q 2 dependence of 
the form factors and the differential decay width for the 
neutral D[ channel. A very similar result is obtained for 
the Dq case. 



As in the previous case, the branching fraction B{D\ — » 
Z?*7r ± ) is again 2/3 in our model because D\ — >• D*tt is 
the only OZI-allowed decay channel. 

Table IVII shows the different helicity contributions to 
the semileptonic width of the reactions B + — > Dil + vi 
and B° — > D^l + vi. The most important contribution is 
given by r^. The ratio Tl/Y = 0.75 gives the probability 
for the final D\ meson to have helicity 0. This result is 
in agreement with the values 0.72 — 0.81 obtained in the 
HQET calculation of Ref. [H|. 

Figure [3] shows the q 2 dependence of the form factors 
and the differential decay width for neutral D\ channel, 
in panels (a) and (6), respectively. Again, a very similar 
result is obtained for the charged case. 

The product of branching fractions are 



B{B+ -> Dll+vi)B{Dl -5- D*-tt+) = 2.57 



B{B° 



Dll + i/i)B{Dl -5- D*-tt+) = 2.57 x 10~ 3 , 
D^l+vi)B{D^ -* D*°tt-) = 2.39 x 10~ 3 , 



(43) 



which in this case compare very well with the latest 
BaBar data Q. 



4- Semileptonic B -> D^lvi decay 

The semileptonic decay is studied by reconstructing 
the decay channel D\ — > D^tt~, using the decay 
chain D* — > D°tt for D* meson and D° — > K~n + 
or D + — > K~ir + ir + for D meson. What is actu- 
ally measured is the product of branching fractions 
B(B+ -> Dfl + ui)B(Df -> D-TT+) and B(B+ -> 



D*°l+vi)B(D* -> D*-n+) 
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FIG. 1. Form factors and differential decay widths for the B + — > Dq°1 + vi decay as a function of q 2 . Very similar results are 
obtained for the B° — > D*^l + ui decay, (a): Form factors predicted by CQM. (b): Differential decay width predicted by CQM. 





FIG. 2. Form factors and differential decay widths for the B + — ► D®l + vi decay as a function of q 2 . Very similar results are 
obtained for the B — > D 1 ~Z + ^( decay, (a): Form factors predicted by CQM. (b): Differential decay width predicted by CQM. 



The first step of this decay involves a semileptonic pro- 
cess which can be calculated using Eq. (j30|) . In Table [VlTl 
we show the different helicity contributions to the total 
width. The main contribution is Yl in both neutral and 
charged D\ channels, providing almost 2/3 of the total 
width. The following one is IV, the rest of the contri- 
butions being negligible. Again our ratio Y^jY = 0.67 
is in agreement with the values 0.63 — 0.64 obtained in 
Ref. UM usin g HQET. 

Figure [4] shows the q 2 dependence in the form factors 
and in the differential decay width, panels (a) and (b) 
respectively, for the B + — > D^l+vi decay. Very similar 
results (not shown) are obtained for the B° —> D*^l + vi 





B+ _> Dfl+vr 


B° -»■ D* 2 -l + ui 


Yu 


0.44 


0.44 


Yu 


2.56 x 10" 8 


2.56 x 10 -8 


Y L 


0.90 


0.91 


Yl 


5.27 x 10" 7 


5.29 x 10 -7 


Y s 


1.54 x 10" 6 


1.55 x 10 -6 


Y 


1.34 


1.35 



TABLE VII. Helicity contributions and total decay widths, in 
units of 10 -15 GeV, for the Do meson. 
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FIG. 3. Form factors and differential decay widths for the B + — >• D®l + vi decay as a function of q 2 . The differences with respect 
B° — > D^[l + vi are negligible, (a): Form factors predicted by CQM. (b): Differential decay width predicted by CQM. 



Branching ratio 



r(D°7r+)/r(L»* 7r + ) 
r(D+7r-)/r(D*+7r-) 

r{D+-K-)/r(D M+ ir- 



Exp. A P 


Microscopic 


1.9 ±1.1 ±0.3 1.80 


1.97 


1.56 ±0.16 1.82 


1.97 


0.62 ± 0.03 ± 0.02 0.65 


0.66 



TABLE VIII. Branching ratios for DJ decays collected by the 
PDG [J] and our theoretical results calculated through the 
two strong decay models. 



case. 



The subsequent strong decays which appear are D 2 — ^ 
D*n~ and D 2 — > Dn~ . In Table I Villi we show the strong 
decay branching ratios obtained with the 3 Po and micro- 
scopic models. They are in good agreement with experi- 
mental data [4j. 

Finally, we obtain the products of branching fractions 
for both decay chains considering that the total width of 
the D 2 meson is the sum of the partial widths of D*tt 
and Dtt channels since these are the only OZI-allowed 
processes 



B(B+ ->■ D* 2 °l+vi)B{D* 2 ° -> D+tt- 



B(B+ -> D?l+vi)B{D* 2 u ->■ D*+tt 



1.44 x 10~ 3 
1.48 x 1(T 3 

0.79 x lfr 3 
0.75 x 10~ 3 

B(B° -> D*-l + vi)B(D*- -> D\-) = \ ]f a * ] n ll 



B(B° -> D* 2 -l + Vl )B{D 2 



* -n_ |0.74xl0- 3 
u n )-lri7ftvin-3 



1.38 x 10" 

0.74 x 10" 
0.70 x 10- 



(44) 



where the first one refers to the calculation using the 3 Pq 



model and the second one comes from the microscopic 
model. These results are in very good agreement with 
BaBar data Q. 



5. Summary of the results 

Final results and their comparisons with the exper- 
imental data are given in Table IIX1 Except for the 
D[(2430), the predictions are in very good agreement 
with the latest experimental measurements, Belle for 
D (2400) and BaBar for D x (2420) and 1^(2460). For 
the D[(2430) there is also a strong disagreement between 
experimental data in the neutral case. 



B. B s semileptonic decays into D** mesons 

The semileptonic decays of B s meson into orbitally ex- 
cited P-wave charmed-strange mesons (D**) provides an 
extra opportunity to get more insight into this system. 

The j q = 1/2 doublet, D* (2318) and D s i(2460), 
shows surprisingly light masses which are below the DK 
and D*K thresholds, respectively. These unexpected 
properties have triggered many theoretical interpreta- 
tions, including four quark states, molecules, and the 
coupling of the qq components with different structures. 
As mentioned before, the D s i(2460) meson has an im- 
portant non-qq contribution. 

We have calculated the semileptonic B s decays as- 
suming that the D** mesons are pure qq systems. For 
the £>* (2318) and D s i(2460), which are below the cor- 
responding D(*>K thresholds, we only quote the weak 
decay branching fractions. Concerning the D s i(2460), 
and as shown in Ref. [15J, the l P\ and 3 Pi probabilities 
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FIG. 4. Form factors and differential decay widths for the B + — > DJ°Z + ^ decay as a function of q 2 . Very similar results are 
obtained for the B° — > D%~l + vi decay, (a): Form factors predicted by CQM. (b): Differential decay width predicted by CQM. 



Belle [1] 

(xicr 3 ) 



BaBar [2, 3] 3 P Mic. 

(xl(T 3 ) (xl(T 3 ) (x!0~ 3 



') 



Do" (2400) 

B(B+ ->■ Dfl+v^BiDf -*■ DJ-K+) 
B(B° -» D^-l+vi)B{D^- -» £>°7r-) 



2.4 ±0.4 ±0.6 2.6 ±0.5 ±0.4 2.15 2.15 

2.0 ±0.7 ±0.5 4.4 ±0.8 ±0.6 1.80 1.80 



Di(2430) 



B(B+ 
B(B° 



n'0;+ 



fi'O 



D*-^ 



d?i + ia)b(d; 

D'-l + v{)B{D'- -^D*°h-) 



< 0.7 
<5 



2.7 ±0.4 ±0.5 
3.1 ±0.7 ±0.5 



1.32 
1.23 



Di(2420) 

B{B+ -> D%l + Ui)B(D% -+ D*--k+) 
B(B° -> D-l + is t )B{D- -> g*y-) 



D2(2460) 

B(B+ 
B(B+ 
B(B+ 



B?l+v,)B(D? 
Dfl+ui)B(Df 
Dl°l + Vl )B{Dl° 



B(B° 
B(B° 
B(B° 



D* 2 -l+u l )B{D* 2 

D^l+^BiD*- 

D* 2 -l+ Vl )B{D* 2 



B 



DID 1 -') 



1.32 
1.23 



4.2 ±0.7 ±0.7 2. 97 ±0.17 ±0.17 2.57 2.57 

5. 4 ±1.9 ±0.9 2. 78 ±0.24 ±0.25 2.39 2.39 



D--JT+) 2.2 ±0.3 ±0.4 1.4±0.2±0.2 W 1.43 1.47 

D*~tt + ) 1.8 ±0.6 ±0.3 0.9±0.2±0.2 W 0.79 0.75 

D M 'n + ) 4.0 ±0.7 ±0.5 2.3 ± 0.2 ± 0.2 2.22 2.22 

>D a -K~) 2.2 ±0.4 ±0.4 l.l±0.2±0.1 w 1.34 1.38 

> D* n~) <3 0.7±0.2±0.1 W 0.74 0.70 

>J) (, V) < 5.2 1.8 ±0.3 ±0.1 2.08 2.08 

0.55 ±0.03 0.62 ± 0.03 ± 0.02 0.65 0.66 



TABLE IX. Most recent experimental measurements reported by Belle and BaBar Collaborations and their comparison with 
our results. The symbol (*) indicates the estimated results from the original data using B D , D {- t) . 



change with the coupling to non-qq degrees of freedom. 
What we do here is to vary these probabilities (includ- 
ing the phase) in order to obtain the limits of the decay 
width in the case of the £) s i(2460) being a pure qq state, 
see Fig. \5\ Assuming that non-qq components will give 
a small contribution to the weak decay, experimental re- 
sults lower than these limits will be an indication of a 



more complex structure for this meson. 

For the decay into Z? sl (2536), our model pre- 
dicts the weak decay branching fraction B(B a s — > 
£> s i(2536)^ + ^) = 4.77 x 1CT 3 and_the strong branch- 
ing fractions B(D sl (2536)~ -¥ D*-R ) = 0.43 (0.47) for 
the 3 Po (microscopic) models. The final result appears 
in Table |Xj It is in good agreement with the existing 
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Experiment 
(xlO -3 ) 


Theory 
(xlO -3 ) 


do (2318) 
B(B° S -> D*. 


o(2318)" 


> + *V> 










4.43 


D s i(2460) 
B{B° S -> D s 


i(2460)" 


> + ^m) 










1.74 - 5.70 


D si (2536) 


i(2536)" 


>+^)S(D s 


i(2536) _ 


->L>* 


"i? ) 


2.4 ± 0.7 [4, 5] 


3 P Mic. 
2.05 2.24 



O s * 2 (2573) 

B(B° -)■ D: 2 (2573)-m + ^)6(I>:2(2573)- 
S(P° -> 7?: 2 (2573)- M + ^)£(-D:2(2573)- 
B{B° S -» j: 2 (2573)~M + ^)g(-P s *2(2573)' 



^Po Mic. 



-if ) 


1.70 1.77 


-if ) 


0.18 0.11 


w-x°) 


1.88 1.88 



TABLE X. Our predictions and their comparison with the available experimental data for semileptonic B s decays into orbitally 
excited charmed-strange mesons. 



> 

o 



° 1.5 




decay part of the reaction, we obtain in our model 



Bp- -+ D-K° 



B(D* S - -> D*~K") = 




0.047 
0.030, 



(45) 



where the first one refers to the calculation using the 3 Po 
model and the second one comes from the microscopic 
model. Our final results can be seen in Table. |Xj 
Besides we predict the ratio 



T(D* s2 -> DK) 



r(D* 2 ->DK) + T(D* s2 ^D*K) 




(46) 



FIG. 5. Decay width for the B° -> D s i(2460)-^ + ^m decay 
as a function of the Pi component probability. The sign 
reflects the relative phase between 1 Pi and 3 Pi components: 
-1 opposite phase and +1 same phase. 



experimental data [4|, which to us is a confirmation of 
our former result in Ref. [15J about the qq nature of this 
state. 

In the case of the _D* 2 (2573) the open strong 
decays are DK and D*K, so the experimen- 
tal measurements must be referred to B{B® — » 
D* 2 (2573)-ii+^)B(D* s2 (2573)- -> L»-ii"_ ) and B(B° a -► 
P* 2 (2573)"^ + ^) S ( £) S *2(2573)- -> D*-R ). 

For the weak branching fraction we get in this case 



S(B° -> # s * 2 (2573) 



> + ^m) 



3.76 x 10 . For the strong 



VI. CONCLUSIONS 

Wc have performed a calculation of the branching frac- 
tions for the semileptonic decays of B and B s mesons 
into final states containing orbitally excited charmed and 
charmed-strange mesons, respectively. 

We worked in the framework of the constituent quark 
model of Ref. [lQ|. The model parameters were fitted to 
the meson spectra in Refs. [10lll4|. Our meson states are 
close to the ones predicted by HQS as expected. 

We have calculated the semileptonic decay rates within 
the helicity formalism of Ref. [22J and following the work 
in Ref. [23]. The strong decay widths have been calcu- 
lated using two models, the 3 Po model and a microscopic 
model based on the quark-antiquark interactions present 
in the CQM model of Ref. [13]. 

From the experimental point of view, Belle and BaBar 
Collaborations provide their most recent measurements 
for the B meson in Refs. [lj and @, Q respectively. For 
the B s meson only the product of branching fractions 
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B(B° S -> L> sl (2536)- i i + ^)i3( J D s i(2536)- -> D—A -0 ) 
has been determined [J] using the experimental data on 
B(b ->• B° S )B{B° S ->• L> s i(2536)-/i + ^)S(L> s i(2536)- -> 
D*~K°) measured by the DO Collaboration Q and the 
PDG's best value for B(6 -» 5°) 0. 

Our results for i? semileptonic decays into Z?q(2400), 
L*i(2420) and I?2(2460) are in good agreement with 
the latest experimental measurements. In the case of 
£^(2430) the prediction lies a factor of 2 below BaBar 
data. Note however the disagreement between BaBar 
and Belle data for the neutral case. 

In the case of B s semileptonic decays, our predic- 
tion for the B(B° S -> D sl (2536)-)U + i/ |i )B(D al (2536)- -> 
D*~ K°) product of branching fractions is in agreement 
with the experimental data. This, together with the 
properties calculated in Ref. [l5[, is to us evidence of 

I 



a dominant qq structure for the D s i(2536) meson. We 
also give predictions for decays into other D** mesons 
which can be useful to test the qq nature of these states. 
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Appendix A: Form factor decomposition of hadronic matrix elements 

Here we give general expressions valid for transitions between a pseudoscalar meson Mi at rest with quark content 
qf 1 qf 2 and a final Mp meson with total angular momentum and parity J p = + ,l + ,2 + , three-momentum — \q\k, 
and quark content 9/'<7/ 2 - The transition changes the antiquark flavor. Following Ref. [23] we evaluate V^(|<f|) and 
A^(\q |) in the CQM through the relations 



V»(\q\) = y/2 mi 2E F (-q) NR { M F , A - \q \k\ J*^(0)|Mj, ) 



NR 



A»(\q\) = ^2mi2E F {~q) NR { M F , A- \q\k\J b A Cfl (0)\Mi,0) NR 



(Al) 



For the different cases under study we will have the following. 



Case CT 



A° m) =V2 m i2E F (-q) [ d* pi (^PWPI))*^ 



I E r E h 
4E fi E fl 



P- 



4ir\p\ 
rh-\q\k-p) p- ( - 



m f i +m/2 



Ef 



A 3 (\q |) = y /2m I 2E F (-q) 



X \Pz 



xp ■ 



m f / +m 



f[ + m /2 



\q\k-p 



n r,) ( \? 



m f ,+m /2 



m h 



m fi + mf, 
q\k-p 



■qk\ 



>< J ^(^r } (w)) mr ] ^- 



E t 



m h 



m f[ + m h 




E f[ E h 
4E f ,E h 



"I 2 \q\k-pi\ 

m.f/ -\-rrif 2 I * I " 



i fi + m h 



\q\k-p 



E f{ E h 



m/ 2 



l f , + m h 



Pz P 




(A2) 



9\-Pi 
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Ef and Ef 1 are shorthand notations for Ef (— — - * — \q\k — p ) and Ef t ( — - ' 2 m — \q\k — p ) respectively and Ef 

E f +mf. 



2. Case CT -¥ V 



Here we have to distinguish two different cases that depend on the total spin 5* of the quark- antiquark system. 



i) Case 5 = 



V, 



(1+,S=0)0, 



A=0 



q\) = - iV3^2 mi 2E F (-q) / d 



E f ,E h 1 



1>XI ±E r E h Airp V^f- 



?(M F (1+,S=0)) 



(P) 



^(Mj(O-)) 
/1/2 



(IP- 



m /2 



m /{ + m h 



V, 



(1+,S=0)1, 



A=-l 



=i\ -y/2mi2E F {-$) / d : 



Sf-^b- 



Itflfcl) p 2 



-Efi-Efi 1 



m f / +m 



y.Z-t-m/j 



!L—\q\k-p) ■ ( mf "+ 2 mf2 W\k-p 



E fi E h 



?(M F (1+,S=0)) 



m/ 2 



m /{ + ™V 2 



±EftE h Airp vKf* 
\q\k\)pl[J- + J- 



(P) 



V, 



(1+,S=0)3, 



A=0 



= -iV3V2m / 2EH-«) / ^ #^A(« (1+ ' S=0)) (P) 



4=-l =0)1 (kl 



^(Mj(O-)) 
/1/2 



(IP- 



m// + m/ 2 



4E fi E fl <±irp vKh 

\q\k\)p 



( "? 2 |g|-Pz "? ; |g|+P, 

mff-\-mf 2 '^ ' -^^ m */ -\-m,f 2 '^ ' -^^ 



V 



£,■ 



%{ 



'/fv^^M)/^^^ (*«'"' 



,S=0)) 



(p) 



;(m 7 (o-)) 
'/1/2 



(IP- 



m/ 2 



i// + m/ 2 



M 






(A3) 



ii) Case 5 = 1 
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v& s ^°(W\ 



E fl E h 1 /?(M F (1+,S=1)) 



=i\ - \/2m I 2E F (-q) / d 3 p\ , ,■..,., 

W 2 V J V AE fi E h 4lT P ^ * 



(p) 



^ ( i + -f =1)1 (kl 



i(Mj(0-)) 
"/1/2 

.>/3 



(|p- 



m/ 2 



m/> + m/ 2 



171*1 



|g|(K-P ) 
E f' 1 E h 



?- v ?2m I 2E F (-q) / rf 



2 

2(M/(0-)) 
/1/2 



E f[ E h 1 /?(M F (1+,S=1)) 



^E fi E h Att P Vhh 



(p) 



V 2 



(|P- 



'/2 



m/j + m/ 2 



171*1 



Pv+^+P»|g| m/ r+L /2 P'+P'-P 2 |7l;^7W 



£ 



/; 



£ f 



^ ,S=1)8 (I7I 



4^r }i ( 



=iJ-^2m,2E F {-q) / rf 



^/(^/i 1 /?(M F (1+,S=1)) 



x«f _)) (l^- '"' 



m /{ + m /2 



' ±E r E h Airp Vfih 

171* I) (p*+p*)( J- -J- 

\ E h E f[ 



(P) 



■ >/3 



=i^-^2mi2E F (-q) 



#pJ E K E * l 



AEf i E fl 4irp 



V Jl -' 2 / J1J2 TOf + mf. 



x < p 2 



1- 



'/{ 



m/ ,+m /2 



171* -P 



m /( + m /2 



|7|*-P 



E fl E h 



rrif> + m h 



m h -m fi p^\g\ 
m f[ +m h E f[ E fl 



(A4) 



3. Case CT -¥ 2 4 



Here we have to distinguish between L = 1 and L = 3. 



i) Case L = 1 



18 



^f =1)1 (k1 



■^ /o OF ? ^ f\l3 E li E h 1 /?(M F (2+,L=1)), 

.. — y /2m I 2E F (- q )Jdp^ w - E -—l<- 



m {2+ ' L ^y 



€ir\\?- 



m h 



m fi + m h 



W\k\) 



'pl-P 2 z-pM ™+ mh P 2 y -Pl+Pz\q\ 



m f[ + m f2 



E 



n 



E h 



4 2 =Y =1)0 (i« 



\)=-^=J2m I 2E F {-q) jd 



E f[ E h 1 



?(M F (2+,L=1)) 



(p) 



2(M/(0-)) 



:Ca ;J dP- T Iglfcl) 

P 2 + p 2 _ 2 2 + 2 |~| ■"/, 
fx fy fz fz\H \ rrifi +mj 2 



AEf^E fl AlTp yfifc 

/pl+pl-2pl-2p z \q\—^ 

'1 



'/2 



V 



E 



fi 



E 



h 



Ati^m) =i^y/2 mi 2E F {-q) I d 3 



I E fl E fl 1 
AE^Efr Awp 



(%£*•" "IP))' 



?(M,(0-)), 



J1J2 ruf + m/ 2 i 



VJ 



m /{ + m /2 



|g|fc-p 



'-'■ 



m /;+ m /2 



|g|fc-p 



E f[ E h 



n 2 Jl\^\ m ^~ m S\ 



^PzpI - pl\q\—_ 



E f[ E h 



A 



(2+,L=l)3, 
A=0 V 



q\) = - iV2y/2m I 2E F (-q) / d 6 p 



E f ,E fl 1 



AE fi E h Amp Vflh 



?(M F (2+,L = 1)) 



(P) 



(M/(0-)), 



*cr ; w 



m h 



m f[ + m h 



m\)\Pz 



V{ 



m f / +m 



/; + m /2 



|g|fc-p 



m /( + m /2 



|(f|fc-p 



E f{ E fl 



+ 



E f ,E h 



2pz{ \q\~P; 



m f > + m h 



m/ 2 



i fi + m/. 



-Itfl-P* + 



(Px+Py) \-Pz + 7T( ~ H 9 



(A5) 



ii) Case i = 3 
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r (2+,i=3)l, 

'A=+l \ 



I 

Vs' 



\-r_:'r ,,L (\q\) =^2 mi 2E F {-q) / d 3 p\ 



J1J2 



(\P- 



m/ 2 



E fj E h 1 

q\k\) 



;(M(2+,L=3)) 



(P) 



1 -(p^2pl-3p z ( ™* \q\-p,)) 
It \ \ y \mr> +mt„ // 



E 



^Pz - 2p p z 



""' 2 -W\-Pz)(pI-p 2 v +p 2 z ) 



i f[ + m h 



1 



E 



P 



/i 



2P 2 V + 3 Pz 



l S[ 



i fi + m h 



\q\-p> 



4 2 ir )o (ki 



5p z y - 2p yPz 
/3 



l /l 



m /{ + m/2 



I«|-^)(p|-^+p|)J 



■£#■% 1 



i^ -V2 mi 2E F (-<l) j d 3 Py ^^^ V ^/ : 



(M(2+,L=3) 



(P) 



P 2 /5p 2 



m/ 2 



P \P 



m f[ + m/ 2 



/;+™/ 2 



?|fc|) 
I? I -P. 



5^ 



-1 



VJ 



Pz+P? , P^+P? 



m f' 1 + m f2 



, E h 

Yi\ +Pz 



Efl 



E h 



E n 



A 



(2+,£=3)3, 



A=0 



-y/2Tm2E F {-q) J d 3 



E fi E h 1 /7(M(2+,L=3)) 



J1J2 



(\P- 



m h 



Mf- 



m f[ + mf 
( 



P ]j 4E fi E h 4irp V^ h 
q\k\) 



\P) 



(pI+pI)(^-i 



V \ 9 



1 



\f[ 



Px+Py- \- mi ,+m f2 n\-P 



V 



A { tl L r m (\q I) = - -j=^2 mi 2E F (-q) j ,l :i p\ 



1 X=+1 



E f[ E h i 



E h E fi 



?(M(2+,L=3)) 




■iE fi E h Airp vKf* 



(P) 



1112 



(\P- 



m h 



2 2 

Pl-Py 



J i 



\m) 



iPz 



3p 2 



5p 2 



m f' 1 + m f2 



\q\-Pz 



m f' 1 + m f2 



m-p* 



pI + pI_pI 

p 2 p 2 p 2 



E h E f[ 



2 2 

Px-Pv 



V{ 



m /;+ m f 2 



l<?l~P; 



V 2 



m/ , +m /2 



I? I ~P; 



£/>£/; 



^ ('|f - 1 



™/2 - m /( 
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2p 2 
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